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ABSTRACT 

The naive implementation of finite width effects in processes involving unstable 
particles can violate gauge invariance. For the example of radiative W produc- 
tion and decay, qq' — > ivy, at tree level, it is demonstrated how gauge invariance 
is restored by including the imaginary part of triangle graphs in addition to re- 
summing the imaginary contributions to the W vacuum polarization. Monte 
Carlo results are presented for the Fermilab Tevatron. 



Radiative W boson production and decay at hadron colliders, i.e. the process qq' — > 
£ ± z/7, is an important testing ground for the Standard Model (SM). At large photon trans- 
verse momentum this process allows the measurement of the WW^j three gauge boson 
coupling [0,0. At small transverse momenta of the photon or when the photon is emitted 
collinearly to the final state charged lepton this process needs to be fully understood when 
trying to extract precise values of the W boson mass and width from the Tevatron data: in 
the CDF detector, for example, photon emission effects in W decays lead to a W mass shift 
of -65 ± 20 MeV (-168 ± 20 MeV) in the W -> ev e (W -> fWy) channel §, and a shift in 
the W width of 70 ± 28 MeV |3J . However, initial state radiation diagrams are not included 
in the Monte Carlo program currently used to simulate these effects ||£|]. 

Thus, a calculation is needed which fully describes initial and final state radiation, which 
implements the effects of the lepton mass, and which incorporates the finite W decay width, 
rV. Finite lepton masses mostly affect photon radiation in the collinear region. Finite 
W width effects, on the other hand, play an important role for large angle soft photon 
bremsstrahlung, where the radiative W decay diagrams and W'y production graphs may 
interfere substantially. In Ref. ||, radiative W decays were considered, including a full 
treatment of lepton mass effects. Initial state radiation was not incorporated into this 
calculation. In Ref. g], we have calculated initial and final state radiation effects. However, 
since we were interested in the measurement of the WW'j coupling, our Monte Carlo program 
did not allow the simulation of collinear final state photon radiation and the implementation 
of finite W width effects was inadequate for soft photon emission. 

Part of the problem can be traced to the difficulty of implementing finite W width ef- 
fects while maintaining electromagnetic gauge invariance 0. The full process qq' — > £ Wy 
proceeds via Feynman graphs with one (initial and final state radiation) and two W propa- 
gators (photon emission from the W). Replacing the W propagator factors l/(q 2 — m^/) by a 
Breit-Wigner form, 1 / (q 2 — + imw^w) , will disturb the gauge cancellations between the 
individual Feynman graphs and thus lead to an amplitude which is not electromagnetically 
gauge invariant when finite lepton mass effects are included. In addition, a constant imagi- 
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nary part in the inverse propagator is ad hoc: it results from fermion loop contributions to 
the W vacuum polarization and the imaginary part should vanish for space-like momentum 
transfers. 

In this Letter we demonstrate how this problem is solved by including the imaginary part 
of WW'j vertex corrections in addition to the resummation of the W vacuum polarization 
contributions. The resulting amplitudes respect electromagnetic gauge invariance, are valid 
for arbitrary lepton masses, combine both initial and final state radiation, exhibit good high 
energy behaviour, and the imaginary parts of the propagators are exactly as derived from 
the W vacuum polarization. Even though we are only considering radiative W production 
and decay here, our method can immediately be generalized to t-channel processes such as 
e + e~ — > e + W~v e — > e + i , e du |7|||. We also present some numerical results of our Monte 
Carlo calculation. 

The general structure is best understood by first considering the lower order process 
(qq' — > i~v to be specific) without photon emission. The corresponding Feynman graphs 
are shown in Fig. |I| Finite width effects are included in a tree level calculation by resumming 
the imaginary part of the W vacuum polarization. When neglecting the fermion masses in 
the loops of Fig. |1] the transverse part of the W vacuum polarization receives an imaginary 
contribution 

while the imaginary part of the longitudinal piece vanishes. In the unitary gauge the W 
propagator is thus given by 

—i ( q >M q u \ i q^q v 



q 2 — m^y + ilmllyy(q 2 ) \ q 2 / — ilmlly^(q 2 ) q 2 



g^-^f(l + ijw)) , (2) 



q 2 - m 2 w + iq 2 -f W V m w / 
where we have used the abbreviation j w = Y w /m w . Note that the W propagator has 
received a q 2 dependent effective width which actually would vanish in the space-like region. 



Since here we are only interested in (virtual) W decays we do not include the corresponding 
step-function factor. 

A gauge invariant expression for the amplitude of the process qq' — > i~v^i is obtained by 
attaching the final state photon in all possible ways to all charged particle propagators in the 
Feynman graphs of Fig. [I]. This includes radiation off the two incoming quark lines, radiation 
off the final state charged lepton, and radiation off the W propagators. In addition, the 
photon must be attached to the charged fermions inside the W vacuum polarization loops, 
leading to the fermion triangle graphs of Fig. |2|. Since we are only keeping the imaginary 
part of Il^(g 2 ), consistency requires including the imaginary part of the triangle graphs only. 
This imaginary part is obtained by cutting the triangle graphs into on-shell intermediate 
states in all possible ways, as shown in the figure. 

The full set of Feynman graphs contributing to qq' — *> l~v^ is depicted in Fig. ^. When 
attaching the photon in all possible ways to either one of the fermion loops or to one of 
the lowest order W propagators in all graphs of Fig. [1|, the remaining sum over W vacuum 
polarization graphs restores the full W propagator of Eq. (0) on either side of the triangle 
graph or the WW'y vertex. Hence, after resummation, one obtains the Feynman graph of 
Fig. |3]c where the full WW^ vertex is given by the sum of the lowest order vertex and the 
imaginary part of the triangle graphs, as defined in Fig. |2|. 

For the momentum flow shown in Fig. |2| the lowest order vertex is given by the familiar 
expression 

- ieTf" = -ie + q^g^ - {q x + fc) V° + (k — q 2 ) a g^) . (3) 

Neglecting the masses of the fermions in the triangle graphs and dropping terms proportional 
to k 11 (which will be contracted with the photon polarization vector s*^ and hence vanish 
in the amplitude) the contributions from the four triangle graphs reduce to an extremely 
simple form. Each fermion doublet /, irrespective of its hypercharge, adds i(g 2 /48tt)Tq to 
the lowest order WW'y vertex IV After summing over all fermion species, the lowest order 
vertex is thus replaced by 
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f*" = r?» ( 1 + £ ifl ] = r?" f 1 + • ^) = rf" (i + llw ) . (4) 

This expression corresponds to the full WWj vertex which is represented by the hatched 
blobs in Figs. ^| and || 

By construction, the resulting amplitude for the process qq' — > £~h"y should be gauge 
invariant. Indeed, gauge invariance of the full amplitude can be traced to the electromagnetic 
Ward identity 

= ( iD w)al (?i) - (<fe) • (5) 

Since 

M**" = - 9?9?) - teV 3 - <&<&)) (i + *rw) , (6) 

and 

/ • j-\ \— 1 / \ / 2 2 , • 2 \ / Qa^/A 2 QaQ/3 /f? \ 

K lD w) a p (?) = (9 - m w + 7^ U Q/3 - I - — , (7) 

this Ward identity is satisfied for the W propagator and WW^ vertex of Eqs. (@) and (f|). 

There are other ways to include finite width effects while restoring electromagnetic gauge 
invariance |7j. One possibility is to multiply the complete tree level amplitude by an overall 
factor F(q 2 ) = (q 2 — m 2 ^)/^ 2 — m 2 ^ + imw^w) for each almost on-shell W propagator ||. 
This procedure implements Breit-Wigner propagators in all resonant graphs but multiplies 
non- resonant contributions by an ad- hoc factor. This is phenomeno logically acceptable if 
the non-resonant contributions are sub-leading. For W'j production, however, the "non- 
resonant" graphs of Fig. ||]a,b are dominant near threshold and multiplying them by the 
factor F(s) decreases e.g. the soft photon emission cross section by 30% or more as compared 
to the correct treatment described above. 

Electromagnetic gauge invariance is also obtained by dropping the triangle graph con- 
tributions and using a naive Breit Wigner propagator ||, 



However, this propagator does not correspond to a resummation of W vacuum polarization 
graphs, the imaginary part does not vanish in the space-like region, and the imaginary part 
of the W mass squared in the longitudinal piece of the propagator is entirely unphysical. On 
theoretical grounds this solution of the gauge-invariance problem must therefore be rejected, 
even though, for the problem at hand, numerical results are very similar to the full treatment 
including triangle graphs. 

The modification of the lowest order WW'y vertex in Eq. (f|) looks like the introduction of 
anomalous couplings gj = k 7 = 1+ijw and one may thus worry that the full amplitude will 
violate unitarity at large center of mass energies \/§. While indeed the vertex is modified, 
this modification is compensated by the effective s-dependent width in the propagator. As 
compared to the expressions with a lowest order propagator, l/(s — m^r), which of course 
has good high energy behaviour, the overall effect is multiplication of the analog of the 
Feynman graph of Fig. |3]c by a factor 

s-m 2 w iT w m w 



s(l + i^w) - mw s-mfy + imwT 



Obviously, G(s) — > 1 as s — > 00 and the high energy behaviour of our finite width amplitude 
is identical to the one of the naive tree level result for W'-f production. In fact, the contri- 
butions from the triangle graphs are crucial to compensate the bad high energy behaviour 
introduced by the g 2 -dependent width in Eq. 

What are the results of our implementation of finite width, fermion triangle, and lepton 
mass effects? In Table I we give the fraction of W ± — > i^v events (£ = e, fx) at the 
Tevatron {jpp collisions at y/s = 1.8 TeV) which contain a photon with a transverse energy 
Et{i) > Eril)- We then compare our results with the expectation for final state radiation 
only (in the narrow W- width approximation), as in Ref. |J. Specifically, we consider the 
phase space region 

p T {£) > 25 GeV , i> T > 25 GeV , \rj(£)\ < 1 , |r/( 7 )| < 3.6 

65 GeV < < 100 GeV , (10) 
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where rrijf denotes the transverse mass of the iv system. The phase space region defined 
by Eq. ( [TOD is very similar to that used by CDF || and D0 [T(J to extract the W mass. 
We then vary the requirements on the minimum photon transverse energy, ££(7), and the 
lepton photon separation, AR(£, 7). QCD corrections and structure function effects largely 
cancel in the cross section ratio considered. Finite lepton masses are seen to hardly influence 
large angle photon radiation whereas the inclusive rate for photon radiation in the electron 
channel is about a factor two larger than the result obtained for muons. Initial state photon 
radiation approximately triples the fraction of W boson events which contain a well isolated 
photon (AR(£, 7) > 0.7). It increases the total number of W — > ev (W — > \xv) events 
which contain a photon by typically ~ 10% (~ 20%). The event fractions presented in 
Table I are insensitive to the maximum absolute lepton pseudorapidity, 77 max , in the range 
= 0.6. . . 1.2. 

Finally, we briefly comment on radiative Z production and decay, qq — > Z7 — > £ + £~'~f 
and qq — > Z — > i^i - ^. Since only diagrams with one Z propagator contribute, the gauge 
invariance problems encountered in the W case when incorporating finite lepton masses do 



not arise here, and the calculation is straightforward. Using the results of Ref. ITT 



we 



have calculated radiative Z production and decay including finite lepton masses at hadron 
colliders. The fraction of Z boson events which contains a photon with Erij) > Eril) is 
about a factor 2 larger than the corresponding number in the W case. For AR(£, 7) > 0.7, 
initial state radiation approximately doubles the number of Z boson events with a photon. 
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TABLES 

TABLE I. Fraction of W — > ev and W — > /iz/ events (in percent) at the Tevatron (pp 



collisions at a/s = 1.8 TeV). We consider events in the phase space region of Eq. fllCf ) 
which contain a photon with transverse energy £t(7) > Et(i) f° r two different minimal 
separations between the charged lepton and the photon, AR(£,j). The results obtained by 
considering final state radiation only (in the narrow W width approximation) are shown in 
brackets. 



W ^ ev W —* \iv 



E° T (l) 


AR(e,j) > 


Ai2(e,7) > 0.7 


A J R( A i, 7 ) > 


AR(n,-/) > 0.7 


0.1 


19.3 (17.5) 


3.63 (1.31) 


11.8 (9.8) 


3.62 (1.32) 


0.3 


14.8 (13.5) 


2.56 (0.95) 


8.76 (7.37) 


2.54 (0.95) 


1 


9.42 (8.63) 


1.45 (0.57) 


5.37 (4.62) 


1.44 (0.56) 


3 


4.58 (4.14) 


0.66 (0.24) 


2.55 (2.14) 


0.65 (0.24) 


10 


0.55 (0.36) 


0.13 (0.01) 


0.32 (0.17) 


0.13 (0.01) 


30 


0.014 (-) 


0.012 (-) 


0.013 (-) 


0.012 (-) 
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FIGURES 

FIG. 1. Feynman graphs for the lowest order process qq' — > £~za The resummation 
of the imaginary part of the W vacuum polarization leads to the Breit-Wigner type W 
propagator of Eq. (Q) which is represented by the shaded blob. 

FIG. 2. Effective WW^ vertex as needed in the tree level calculation of W'j production. 
In addition to the lowest order vertex the imaginary parts of the fermion triangles must be 
included (see Eq. ([|)). 

FIG. 3. Feynman graphs for the process qq' — *> i~v^(. The shaded and hatched blobs 
indicate the full W propagators as defined in Fig. [I] and the full WW'j vertex as given in 
Fig. 0, respectively. 
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